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Abstract. Continuous time models in the theory of real options give explicit 
formulas for optimal exercise strategies when options are simple and the price of 
an underlying asset follows a geometric Brownian motion. This paper suggests 
a general, computationally simple approach to real options in discrete time. 
Explicit formulas are derived even for embedded options. Discrete time processes 
reflect the scarcity of observations in the data, and may account for fat tails and 
skewness of probability distributions of commodity prices. The method of the 
paper is based on the use of the expected present value operators. 
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1 Introduction 

The goal of the paper is to provide a general framework for pricing of real options 
in discrete time; since option pricing cannot be separated from the problem of 
the optimal exercise time of an option, we solve for the optimal timing as well. 
Valuation of real options and optimal exercise strategies explicitly derived in 
the paper are relevant to many practical situations where individuals make at 
least partially irreversible decisions. Simple real options cover such situations as 
timing an investment of a fixed size, the scrapping of a production unit, a capital 
expansion program, etc. Embedded real options are relevant in such cases as new 
technology adoption, timing an investment (partially) financed by debt with an 
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embedded option to default in the future, human capital acquisition, and many 
others (see, for example, papers in the volume [Sj, and the bibliography 
therein). Usually, continuous time models are used. We believe that the discrete 
time approach has certain advantages over the continuous time approach to 
modeling of real options. 

The discrete time approach allows one to incorporate in a more tractable way 
such important aspects of economic reality as "time-to-build" . Some options - 
for example, renewing a labor contract or awarding a grant to a research project 
- can be exercised only at certain dates (typically, once a year) . Further, the data 
in economics is much more scarce than in finance, where continuous time models 
were originated and developed. Observations of some variables are available only 
quarterly or even annually, and it is easier to fit a transition density to the data 
in a discrete time model than in a continuous time one, when only a list of 
observations at equally spaced moments in time is available. There are also 
natural questions, such as the rates of job creation and job destruction (see, e.g., 
[T^ and JO]), which cannot be addressed in the Brownian motion setting under 
the standard assumption that labor is instantly adjustable, because the rates will 
be infinite. In the discrete time model, one obtains finite closed form solutions. 
Finally, discrete time modeling is natural for economics, where agents do not 
make, say, investment decisions every instant unlike the investors in financial 
markets. 

In fact, the appearance of continuous time models in economics was mainly 
due to the success of these models in finance, and their tractability. Notice, 
however, that even in finance, the success of simple models is gone. Relatively 
simple models use the (geometric) Brownian motion as the underlying stochastic 
process (see, for example, |E1)- However, these models proved to be rather inac- 
curate, and have been amended in many ways, none of which being as tractable 
as the Brownian motion model. In particular, nowadays processes with jumps 
are widely used (see, for instance, and ^Hl)- The normality of commodity 
price processes is rejected by the data as well (see, for example, J3] or ^[). 
Assume, for the moment, that we are willing to sacrifice the accuracy of a model 
for its simplicity, and decide to use the Brownian motion model. The apparent 
advantage is a well-known manageable scheme: with the help of Ito's lemma 
write down a second order differential equation for the value of an option, em- 
ploy economic arguments to add appropriate boundary conditions, such as value 
matching and smooth pasting, and, using the general solution to the differential 
equation, reformulate the problem as a system of algebraic equations. In ele- 
mentary situations, a closed form solution can be derived, and a simple exercise 
strategy results; in other situations, numerical procedures are available. 

Unfortunately, this scheme is uncomplicated only when one considers really 
simple options. If a part of the option value comes from instantaneous payoffs 
due at certain future dates, a closed form solution to the optimal exercise problem 
is no longer available. If we consider an option whose value comes from different 
streams in different regions of the state space, the resulting system of algebraic 
equations may involve too many unknowns, and it becomes messy indeed. If we 
consider embedded options, then there is no general result about the optimal 
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exercise rule, apart from the heuristic smooth pasting condition, and it is not 
even clear whether the formal solution satisfying the smooth pasting condition 
exists. In addition, if one incorporates jumps, the intuitive justification for the 
smooth pasting condition in is lost, and there is no reason to believe that 
this principle always holds^. In fact, it may fail as it was shown in [S] and 

From our point of view, the following two drawbacks of the standard ap- 
proach are conceptually even more important. First, the formulas provided by 
the standard approach are certain analytic expressions, which have no clear 
economic meaning by themselves, and the optimal exercise boundaries are just 
numbers. The classical example is the investment rule formulated in |14j in 
terms of the correction factor to the Marshallian law: the correction factor is 
an expression in terms of a root of a certain quadratic polynomial (fundamen- 
tal quadratic) . Second, is it really necessary to write down the equation for the 
value function, the boundary conditions, and solve the (free) boundary problem? 
Indeed, the option value by definition is the expected present value (EPV) of an 
instantaneous payoff or stream of payoffs. Suppose, we have formulas for these 
expected present values for model payoffs. Then using the argument similar to 
the backward induction in discrete time models with the finite time horizon, we 
can calculate the option value step by step. We solve for the value and optimal 
exercise time (if necessary) of the most distant option in a sequence of embedded 
options, then we move to the second to last option, etc. 

This paper provides a general framework of this sort, with explicit formulas 
for each model situation, under the assumption that the underlying stochastic 
process is a random walk (under a risk-neutral measure chosen by the market), 
X(, on the real line, that is, Xt — Xq + Yi + Y2 + ■■■ + Yt, where 11, 12, . . ■ are 
independently identically distributed random variables on a probability space O, 
and Xf) is independent of Yi, I2, • • ■ This process specification implies that the 
dates when observations and/or decisions to exercise options can be made are 
equally spaced, and time periods are normalized to one. We calculate the EPV's 
in the following model situations: (1) an instantaneous payoff tomorrow; (2) a 
perpetual stream of payoffs; (3) an instantaneous payoff which is due when a 
barrier is crossed from below (or from above) ; (4) an entitlement for a stream of 
payoffs which comes into effect when a barrier is crossed from below (or above); 
(5) a stream of payoffs which will be lost when a barrier is crossed from below 
(or from above). Using these EPV's, we solve for the optimal exercise time of 
an option, calculate the expected waiting time till an option will be exercised, 
and find the optimal gradual capital expansion strategy. 

Let g{Xt) be the payoff function, and q G (0, 1) the discount factor. The 
first two cases are simple: the calculation of qE^[g{Xi)], and the summation of 
an infinite series X^t^o I'^-^^idi-^t)], respectively. For g{Xt) a linear combination 
of exponents, the calculations in the first two cases reduce to the calculation of 
values of the moment generating function of Yi, and closed- form solutions result. 
These observations are fairly standard and commonly used. The novelty of our 
approach is that the solutions in all the cases above are expressed in terms of the 

^We are grateful to Avinash Dixit for pointing out to us that this principle may not hold 
in a discrete time model. 
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expected present value operators (EPV-operators) . For a random walk X, and 
a given stream of payoffs g{Xt), the EPV-operator (denoted by U\) calculates 
the EPV of the stream g{Xt): 



= E 



t=o 



q'giXt) \Xo^x 



This operator gives the EPV in case (2). Cases (3) and (5) reduce to cases (2) and 
(4). To express the EPV of a payoff in case (4), in addition to the EPV-operator 
of X, we need the EPV-operators of the supremum process Xt — maxo<s<t Xg 
and the infimum process -X_ = mino<s<( X^. These EPV-operators act as follows 
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Let l[/i,+oo) denote the indicator function of the interval [h, +oo) and the mul- 
tiplication operator by the same function. We prove that the EPV of a stream 
g{Xt) that accrues after the barrier h is crossed from below is the composition 
of the EPV-operators U% and U^, and Ir^, 



V{h;x) ^ {1 - q)Ul^l^h,+oo)U],g{x) 



(1.1) 



Similar formulas are obtained for the other expected present values listed in 
cases (3)-(5). 

Assume that both the supremum and infimum processes are non-trivial, i.e., 
jumps in both directions happen with non-zero probabilities. Then the form 
of the solution (|1.1|) and simple general properties of the EPV-operators allow 
us to give a very short (approximately a page) proof of the optimal exercise 
rule: exercise an option which yields a stream of increasing payoffs g{Xt) the 
first time when the EPV of the stream g under the infimum process becomes 
non- negative . In this case, the record setting low payoff process matters. This 
reflects Bernanke's bad news principle spelled out in T. Similarly, for a decreas- 
ing payoff function g, the record setting high payoff process matters: exercise the 
option the first time when the EPV of the stream g under the supremum process 
becomes non-negative. Notice that a sufficient condition for optimality makes 
a perfect economic sense: the payoff function g{Xt) is a monotone function of 
the stochastic factor Xf. Record setting news principles were obtained in jj] 
in the continuous time model, for a special case of payoff streams of the form 
Ge^ — C(l — q) or C(l — q) — Ge^, by using the well-known results for the per- 
petual American call and put options, respectively. The present paper extends 
the record setting news principles to general continuous monotone streams in 

•^For instance, in timing an investment of a fixed size, g is the stream of revenues net of 
(1 — q)C, where C is the fixed cost of investment. 
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discrete time models. If the payoff stream, g, is an increasing function of the 
stochastic factor Xt, then it is optimal to exercise the right for (respectively, to 
give up) the stream when the EP V of the stream calculated under the assumption 
that the process {Xt} is replaced by the infimum process becomes non- negative 
(respectively, non-positive). For a decreasing payoff function g{Xt), one only 
needs to substitute the supremum process for the infimum process in the above 
statement. 

The reader may think that the calculation of the action of the expected 
present value operators is difficult. For a general random walk, this is really the 
case. However, one of the advantages of the discrete time setting is that the 
transition density of a random walk can be approximated by exponential poly- 
nomials on each half-axis with desired accuracy and simplicity, although there 
is certainly a trade-off between the two. The family of transition densities given 
by exponential polynomials is fairy flexible, and such densities can account for 
fat tails and skewness observed in empirical distributions of commodity prices. 
For a transition density of this family, and g{x) given by exponential polynomi- 
als (different exponential polynomials on different intervals of the real line are 
admissible), all the functions Ujrg{x), U'^g{x) and U'^g{x) are of the same class 
as g{x). Therefore, we can use general formulas derived for model situations 
several times (embedded options), and each time we will calculate the integrals 
of a simple structure. All the tools needed for the calculation of these integrals 
are the fundamental theorem of calculus and integration by parts. It is unnec- 
essary to write down differential equations and solve them. In fact, after the 
roots of a certain polynomial are found (we will call it the fundamental ■polyno- 
mial by analogy with the fundamental quadratic in the Brownian motion case: 
see mi), all the calculations reduce to straightforward algebraic manipulations. 
Even in the case of embedded options of complicated structure, there is no need 
to consider systems of algebraic equations, and timing an option reduces to the 
calculation of a (unique) zero of a monotone function. At the same time, the 
method of the paper allows one to solve the optimal exercise problems for much 
wider classes of payoffs than the standard approach does. 

The rest of the paper is organized as follows. In Section 2, we present general 
formulas for EPV's in model situations (l)-(5) and determine the optimal exer- 
cise time and the waiting time till an option is exercised. In Section 3, explicit 
formulas for the results of Section 2 are presented. In Section 4, timing a capital 
expansion program is solved by reduction to a sequence of investment opportu- 
nities of fixed size, and formulas for the optimal capital stock and the value of 
the firm are derived. The last formula shows clearly that the value of the firm 
does not evolve as a Markovian process. A similar result can be obtained in the 
continuous time model. Therefore the value of a firm cannot be modeled as a 
Gaussian process as it is often done in finance. Technical details and proofs are 
given in the appendices. 
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2 Model situations 



2.1 Instantaneous payoff and payoff stream 

The standing assumption of the paper is that the underlying stochastic process 
Xt is a random walk on R, which has the transition density, p (the method of 
the paper can be applied to random walks on lattices as well). The transition 
operator, T, is defined by 

(T/)(x) = E-[f{X,)] = E[f{X,)\ Xa=x]. 

Given p, one calculates the EPV of a stochastic payoff tomorrow: 

/+CX; 
p{y)f{x + y)dy. 
-oo 

To compute the EPV of a stochastic payoff t periods from now, we use the 
Markov property of a random walk: 

E^[q'g{Xt)]^q\T'g){x). 

The next step is to calculate the EPV of a stream of payoffs: 



U],g{x) = E^ 



^Y.q\T'g)ix). 



Under the condition q G (0, 1), the operator / — gT is invertible in L°°(R), and 
the inverse is (/ — qT)^^ ~ X^t^o^*^*' hence 

Uj,g{x)^{I-qT)-^g{x), (2.1) 

for a bounded measurable g. Under additional conditions on the transition 
density, 1)2.1(1 can be extended to unbounded g. From 12.1|l . we conclude that 

U%{I-qT) = {I-qT)Ul=I. (2.2) 

This equation allows one to express an instantaneous payoff G{Xt) in terms of 
the EPV of a stream and vice versa: 

g{Xt) = {I-qT)G{Xt), (2.3) 
G{Xt) - Ulg{Xt). (2.4) 

Clearly, equation H2.4|) is the value, at date of the stream of payoffs that starts 
at date t as well. To evaluate such a stream at the initial date, we compute 

E^[q'G{Xt)]=qVU'^g{x). (2.5) 



6 



2.2 Payoff stream that starts to accrue at a random time 

The list of the EPV's presented in the previous subsection conveys no new infor- 
mation to any person famihar with the basics of economics or finance. Suppose 
that we want to price the stream of payoffs that starts to accrue after the under- 
lying stochastic variable, Xt, crosses a certain barrier, h. In this case, we need 
to replace the deterministic time, on the LHS in (|2.5() with a random time, 
T, and (|2.5|) is no longer valid. Fortunately, it is still possible to compute the 
value E^[q'^ /{Xt)] in terms of the EPV's of some payoff streams. To do this, we 
need to distinguish between two cases: Xt crosses h from below, and Xt crosses 
h from above. Denote t"^ = Th — niin{t \ Xt > h}, and consider 



V+ix;h) =E'' 



E I'aiXt) 



.t = T + 



Under natural conditions on g, 



V+ix; h) = E'- U%g{XT^) = (1 - q)U''^l[h^+^)U%g{x) 



(2.6) 



Equation H2.6|l holds due to the quite deep result, known as the Wiener-Hopf fac- 
torization formula, that establishes the relationship among the EPV-operators 



[/|., [/|, and [4: 



Ul = (1 - q)UWl 



'X — y-^ ^ 'ij'^x'^2L 

(for details, see Subsection 3.2). The Wiener-Hopf factorization formula (|2.7|l 
per se has no economic intuition behind it (and neither does Ito's lemma, both 
being general deep mathematical results), but using this formula, one obtains 
values of real options and optimal exercise strategies in economically meaningful 
terms. 

Similarly, for the case of crossing from above, set — — min{t \ Xt < h} 
and 



V-{x;h) = E"" 



OO 

E 



q*9iXt 



Then the value of such a stream is given by 



V-{x;h) = E- [g^"C/^5(X,-)J = (l - q)U'^l^_^^h]Uy{x). (2.8) 

It is easy to see that H2.8|l is just a version of H2.6|l : change the direction on 
the real line, then the supremum process becomes the infimum process, and vice 
versa. 

Notice that H2.6|l and H2.8() price the streams of payoffs that are described in 
case (4). The model situations in case (3) reduce to these formulas if we express 
an instantaneous payoff as the EPV of a stream using H2.3|l : 



E"" 



cT'CiX,-] 



= [l - q)Umn.+oo)UW - qT)G{x); 
= (l-g)Ull(_oo,/^]^7|(/-gT)G(a;). 
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2.3 Payoff stream that is lost at a random time 

The case when the right to the stream g{Xt) is lost when Xf crosses h from 
below is reduced to the cases considered above: 



V+(x;h) 



t=0 



E ^'aiXt) 



t=T + 

= Ulg{x) - (1 - q)U'il[n,+oo)Ulg{x). 
The Wiener-Hopf factorization formula (|2.7I) allows us to proceed: 



\x-h) = {l-q)Ul(ulg{x)~^n.+oo)U'k9{x)) 
= {l-q)Um^^^n)Ulg{x). 



(2.9) 



In a similar way, if the right for the stream is lost on crossing h from above, 
we derive 



V-(x]h) = E'' 



E ^'aiXt) 



t=o 



= (1 - q)Uxi(h.,+oo)U]^gix). 



(2.10) 



2.4 Optimal exercise time 
2.4.1 Time to enter 

Suppose now that g is the stream of payoffs specified for a real option, for 
example, an option to invest capital C into a technology. If an entrepreneur 
chooses to invest at time r, then the technology will produce a commodity at 
rate G, starting from time t + 1, ever afterward. Investment is irreversible. One 
may view C as the present value of a deterministic stream of expenditures, to 
which the investor commits at time r. By (|2.3|l . this stream is c = (1 — q)C. 
The output is sold on the spot at the market price Pt = e^* . Obviously, we are 
facing the situation to which p. 6(1 apphes with g{Xt) = qGEt[e^*+'^] — (1 — q)C. 
In Appendix A, we show that if g is non-decreasing and U'^g{x) changes sign 
only once, the optimal exercise (log) price (the optimal investment threshold in 
our example), h* is a unique solution to the equation 



wix) = iUj,g){x) = 



(2.11) 



We stress that contrary to the Brownian motion model and other models with 
continuous trajectories, the exercise boundary in the present model is not nec- 
essarily the price at which the option is exercised. It is optimal to exercise the 
option the first time r such that Xr > h* . For the case of irreversible investment, 
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H2.11|l is equivalent to 



(gGe^'+^ - (1 - 



,t=o 



0. 



(2.12) 



The last equation says that investment becomes optimal when the EPV of the 
project calculated under the assumption that the original price process is re- 
placed by the infimum process becomes non-negative. If the investment thresh- 
old is chosen optimally, then the value of the option to invest, as given by H2.5|l . 
is 

V+{x; h*) = (1 - q)U]^l[h^,+oc)Ulg{x). 
We use (|2.12ll to write 



C^e" E 



,t=i 



Therefore, 



Y^q' (<7Ge^*+i - (1 - g)C) 



J2 q*Ge^' I ^0 = 



,t=i 



(2.13) 



C 



and 



V+ix;h*) 



= {e - e'* )Ce- 
-'^*(l-g)([/|l[,.^+oo)(-)(e--e''*))(x) 



Ce 

Ce-^'{\-q)E 



where (e"^' — e'' )+ = maxKe'^* — e'* ),0}. Now we can express Tobin's Q as 



Q{x) 



y+(a;;/t*) 
C 



E^ 



ei^'/il-q) 



Hence Tobin's Q is the ratio of the EPV of the stream of payoffs {e^* — e'' )+ 
and the EPV of the perpetual stream e'* . Notice that the investment thresh- 
old is determined by the infimum process, and Tobin's Q is determined by the 
threshold and supremum process. 

Alternatively, using (|2.13|) . we can write the value V^{x; h*) as 



V+{x;h*) = {l-q)E 



E<?*Ge^* |Xo = 



The last representation factors out contributions of the infimum and supremum 
price processes to the value of the option to invest. The first expectation on the 
RHS decreases as the probability of downward jumps in prices increases, and 
the second expectation increases with the probability of upward jumps. If both 
probabilities increase, the overall effect on the value V'^{x; h* ) is ambiguous. 
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2.4.2 Time to exit 



The option considered above is similar to the perpetual American call option. 
In this subsection, we consider the optimal exercise strategy for a generalization 
of the perpetual American put option, which is relevant to irreversible decisions 
such as scrapping a production unit or exit; it also applies to the optimal timing 
of default on a debt. 

Consider a firm in a deteriorating environment. The firm's profit, f{Xt), is 
falling, on average, so that it may become optimal to discontinue operations at 
some point in time, and sell the firm's inventory for the scrap value C. If the 
firm's manager makes a decision to scrap the inventory at time r, the operations 
will be discontinued at time r + 1, and the scrap value will be available at time 
T + 2. When the firm discontinues its operations, it looses the stream of profits 
whose EPV is 



.t = T + l 



The scrap value, C, can be viewed as a deterministic stream of payoffs, c = 
(1 — g)C, that starts to accrue at r + 2. Hence, at time r, the option value 
of exit is U'^g{Xr), where g{x) = q^{l ~ q)C — qTf{x). If the decision to exit 
is made when the underlying stochastic factor, Xt, crosses a barrier h from 
above, then the option value to exit is given by (|2.8|) . Suppose that g{x) is 
non-increasing and U'j^g{x) changes sign only once. Then the optimal exercise 
price of the option is a unique solution, call it /i* , of the equation 

w{x)^{Uy){x)^0 (2.14) 

(see Appendix A). For g in the case under consideration, (|2.14|) is equivalent to 



Y,q\q^{l-q)C-qf{Xt+,)) 



.t=0 



(2.15) 



or 



f{x)+q'C^U]^f{x). (2.16) 

The LHS in H2.16|l is the EPV of the instantaneous payoff received if the opera- 
tions are terminated, the RHS is the EPV of the stream of profits evaluated under 
the assumption that the original stochastic process is replaced by the supremum 
process. Rule (|2.16|) says that it is optimal to exit when the difference between 
the former and the latter EPV's becomes non-negative. 
Suppose that f{x) — Ge^, then from 12.15II . one obtains 



q^C ^ e'^'E 



whence 



Y,q\q\l-q)C~qf{Xt+,)) 



.t=0 



= q2^e-''*(e''- -e"). 
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and from (|2.8|l . we derive 

V-{x;K) = 9'Ce-''-(l-9)(c4l(-oo,/..](-)(e''* -e)) (x) 



q^Ce-'''{l-q)E'' 



J2 q'Ge^' I Xo = 



Lt=0 



2.5 Expected waiting time 

We consider timing the investment of a fixed size once again. Assume that 
the spot log- price x is less than h* , and consider the waiting time Rx till the 
investment is made. This is the random variable defined by 

Rx = min{< >Q\Xt> h*}. 

The expected waiting time can be calculated as follows: 



E\Rx 



'^(-oo,h*){xt) 

.t=0 



hm E"" 



.t=0 



lim iE'' 



.t=0 



lim E"" 



,t=0 



= ^Hm^{(l-g)-i-t/|l[,.,+^)(x)} 



(2.17) 



Notice that here one must consider the process under the historical measure 
and not a risk-neutral one. Explicit formulas for the expected waiting time are 
presented in Appendix C. 

3 Explicit formulas for model situations 

In this section, we show how to calculate the EPV when the transition density 
is given by exponential polynomials. 



3.1 Instantaneous payoff and payoff stream 

3.1.1 Calculation of the expected present value of an instantaneous 
payoff at a deterministic moment in the future 

Let p be the transition density of the random walk Xt = Xq + Yi + --- + Yt, and 
let 



M{z) ^ E[e 



zYi] _ 



+ 00 



e'"'p(x)dx 



be the moment generating function of Yi. If G is an exponential function: G{x) = 
e^^, then the transition operator, T, acts as a multiplication operator by the 
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number M{z), and we obtain qE^le^-^^] — qTG{x) = qM{z)e^^. Similarly, one 
easily calculates the expected present value of a payoff G{Xt) at a deterministic 
moment t in the future: g*i;^[e^-^*] = {{qTfG){x) = (gM(z))*e^"=. Equations 
\1.'6\ and (|2.4|) can be used to represent a payoff G as the expected present 
value of the stream g{Xt) = ((/ - qT)G){Xt); G{Xt) = W^giXt). In a special 
case of an exponential function G{x) = e^^, with z satisfying the condition 
1 - qM{z) > 0, we have g{x) = (1 - qM{z))G{x). 



3.1.2 Expected present value of an ongoing project: capital and labor 
are fixed 

Consider a firm that produces a commodity which is sold on the spot at the 
price Pt = e^* . We impose the following restriction on the price process: 

qM{l) < 1. (3.1) 

Restriction H3.1|) is appropriate when capital, K, and labor, L, employed by the 
firm are fixed or capped. Indeed, for K and L fixed, the EPV of the revenue 
stream grows each period by factor M(l) = E [e^^], and it is discounted back 
at rate q. Hence, the EPV of the stream of revenues is given by 

PoG(X,i)f:(,M(l))* = ^^^§^, (3.2) 

where G{K, L) is the production function of the firm. For the series in ijX^ to 
converge, it is necessary and sufficient that 1)3. 1(1 holds. 

3.1.3 The case of costlessly adjustable labor 

In other situations, more stringent conditions than (|3.1|l may be necessary. For 
instance, consider a firm with the Cobb-Douglas production function G{K, L) = 
dK^L^~^, which faces the fixed labor cost w, and can instantly and costlessly 
adjust labor. A similar situation was considered in W for a continuous time 
model of irreversible investment. At a given price level of the output, P, the 
firm chooses L to maximize PdK^L^^^ — wL. From the F.O.C. 

(1 - e)PdK^L-'^ - w = 0, 
we find L = K{{1 — 9)Pd/wY^^ , and hence the revenue is i? = AP^/^ , where 

We conclude that the EPV of the revenue stream 

°° AP 

APoY.i,Mimr-y^^^ (3.3) 
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is finite iff qM(l/0) < 1. Notice that if increases in prices are anticipated, that 
is, M(l) > 1, then M{l/9) > M(l) (recall that the moment generating function 
is convex), and therefore, the condition qM{l/9) < 1 is more stringent than 

3.1.4 The case of stochastic prices and operational costs 

One can allow for a stochastic operational cost as well, and the model remains 
quite tractable provided both prices and variable costs depend on the same 
stochastic factor, call it Xt- A simple example (with fixed capital and labor) 
would be R{Xt) = e^'G{K,L), and C{Xt) = a + be^^^ where a, 6, 7 > 0, so 
that the profit is 

g{Xt) = R(Xt) - C{Xt) = e^'G(if,L) ~a~ 6e^^*. 

These examples show that it may be necessary to consider payoffs of various 
kinds; hence, some general regularity conditions on g are needed. 

3.1.5 Calculation of U'j^g{X) for general payoffs 

The EPV of a stream g{Xt) can be easily calculated when g is an exponen- 
tial polynomial, that is, a sum of products of exponents and polynomials. For 

g{Xt) = e^^', the result is 

f/^e^^ = E- 

provided 1 — qM{z) > 0. By linearity, H3.4|l extends to 5 a linear combination of 
exponential functions, and more generally, for integrals of exponential functions. 
Hence, if an explicit formula for the moment generating function is available, 
U\g{X) can be calculated for all payoff streams g of interest. 

Sufficient conditions for U'^g{x) of a measurable stream g{Xt) to be finite 

are 

\g{Xt)\ < Ccxpia+Xt), > 0, (3.5) 
\g{Xt)\ < Cexp{a-Xt), Xt<0, (3.6) 

where constant C is independent of Xt, and satisfy 

l~qM{a^)>0. (3.7) 

These conditions are necessary if g is monotone on each half-axis. 

If we confine ourselves to a wide and fairly fiexible class of probability den- 
sities given by exponential polynomials on each half-axis, then for a general g, 
the calculation of U'^g{X) reduces to simple integration procedures. Consider 
first the case when the transition density is of the form 

p{x) = c+X+e-^^n[o,+oo){x) + c-i-X-)e-^'n^_^^o]{^), (3.8) 



t=o 



1 -qM{z)' 



(3.4) 
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where c^,c > 0, and A < < A"*". If we want to have a continuous p, we 
must require that c+A+ + c^A^ = 0, and then the normahzation requirement 
Af(0) = lleadstoc+ A-/(A" - A+), = A+/(A+-A-). We have constructed 
a two-parameter family of probabihty densities. The moment generating function 
is 

1 1 



C+A+ c-\- -A-A+ 

M{z) = — + — = — 

\+ — z A — z A+ — A 



A+ - z A- 



It is easily seen that 1 — qM{z) has two real roots: /3 G (A ,0), and (3'^ G 
(0, A^), which are the roots of the quadratic equation 

z2 - (A+ + A-)z + (1 -9)A+A- = 0. (3.9) 

We find 

/3± = 0.5 • (a+ + a- ± v/(A+ + A-)2-4(l-9)A-A+) , (3.10) 

and represent 1/(1 — qM{z)) in the form 

1 a+ 



1 - qM{z) (3+ ~ z 13- 



(3.11) 



where a± = {(3^ ~ z)/{l ~ qM{z))\^^p± = l/qM'{(3^). Set g{x) = e^^. Using 
(|3.8|l and 13.11|l . it is straightforward to check that the RHS in (|3.4I) is equal to 
the RHS in the equation below 

/•oo 

U%g{x)=a+ e-^*yg(x + y)dy-a- e'^^ygix + y)dy. (3.12) 



oo 



By linearity, H3.12|l holds for linear combinations of exponents, and more gener- 
ally, for wide classes of functions which can be represented as integrals of expo- 
nents e^^ w.r.t. z. One can also use more than one exponential on each axis, 
and obtain more elaborate probability densities (see Appendix B for details). 

3.2 Payoff stream that starts to accrue at a random time 

As an example, consider the problem of timing an investment of a fixed size (see 
Subsection 2.4.1). We need to compute the value given by (|2.t)f) . where 

g{Xt) = qGEt[e''*+'] - (1 - q)C = gAf (l)Ge^' - (1 - q)C. 

The operator {l — q)U'^ acts on an exponential function e^^ as the multiplication 
operator by the number (1 — q)/{l — qM{z)): 

1-9 ^zx 



[l^qWle^ 



^ l-qM{zy 

Similarly, the (1 — q)U'^ and (1 — q)U'^ act on an exponential function e^^ as 
multiplication operators by certain numbers, call them k^(z) and k^(z): 

(l-g)^7|.e^- = <(z)e^^ {1 - qW^e/^ = {z)e^\ (3.13) 
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These numbers are 



and 



«-(z) = (l-g) (Ule/ 



) =(1-9)^ 



J2 q'gi^t) I Xo = 



>*.9(^t) I ^0-0 



respectively. The Wiener- Hopf factorization formula (see, e.g., JH]) states that 

1-9 



l-qM{z) 



(3.14) 



It follows that (|2.7|l holds: (1 — <z)^^^x- other words, we can factorize 

the EPV-operator U'^ in terms of the EPV-operators U'^ and U'^. Indeed, if 
we apply both sides of (|2.7|l to e^^, we obtain identity on the strength of H3.14|l . 
By linearity, (|2.7(l extends to linear combinations of exponential functions, and 
more generally, to wide classes of functions which can be represented as integrals 
of exponentials w.r.t. parameter z. 

Under an additional very weak regularity condition on the density p (piece- 
wise continuous p are allowed), it is proved in 7 that if g is measurable and 
satisfies the growth condition 13.5|l . then for r+ — Tj^ , (|2.()|) holds. It is possible 
to prove (|2.t)|l without additional conditions on p. 

If the probability density p is given by exponential polynomials on the positive 
and negative half-axis, then all the functions in H3.14I) are rational functions. 
The calculation of the factors k^{z) reduces to the calculation of roots of the 
numerator of the rational function 1 — qM{z), that is, of a polynomial. To 
obtain k+(z), one collects all the factors in the numerator and denominator of 
the rational function (1 — q)/(l — qM{z)), which do not vanish in the half-plane 
5Rz > (the real part of z is non-negative), and k~{z) contains the factors which 
do not vanish in the half-plane SRz < (see, e.g., ^T] and [3- In many cases 
(including the examples above), all the roots are real. Now the same argument 
as at the end of Section 2 shows that the action of operators and is given 
by simple integral operators. As the result, we obtain an effective procedure for 
the calculation of U'^g{x) and U'^g{x). In this section, we consider the case of 
the probability density 13.8|l : for the general case, see Appendix B. 

First, we calculate the the moment-generating function M{z), find the roots 
of the equation 1 — qM{z) — (see (|3.1U|l '). and represent k'^{z) and k~{z) as 
the sums of constants and simple fractions: 



and 



K-{z) 



{X+-z)l3+ _ 13+ /3+(A+-/3+) 

A+(/3+ -z) A+ X+{p+ ~z) ' 

(A- -z)r ^ r r(A- -r) 

X+ip- -z) X- X-ip- -z) ' 



(3.15) 
(3.16) 
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Using H3.13II . 1)3. 15|) . H3.16|l and the same considerations as in the derivation of 
H3.12|l . we obtain for a continuous g satisfying (|3.5|l 



/3+(A+- /3+) 
A+ 



-^^yg{x + y)dy, (3.17) 



and for a continuous g satisfying (|3.()ll . we calculate 
(1 - Q)Ux_g(x) = —9{x) 



e-P'yg{x + y)dy. (3.18) 



Under condition 1)3. 7|l . the roots (3^ are outside [cr~,o'+], therefore 1)3.5(1 and 
H3.6|l ensure the convergence of the integrals in ((3.17() and l(3.18|l . 

Now we return to the investment problem of Subsection 2.4. From (|3.13|l . 
wc find 

Uyix) = (1 - q)-'K-{l)qM{l)Ge- - C. 

Set 

u{x) = l[h^+^){x)iUy)ix) = l[,,+oo)(x)[(l - q)-'K-{l)qM{l)Ge^ - C], 

and calculate the value V'^{x; h) — E'^lq"^^ U'^g{X^+)] of the investment project 
at time (provided h is chosen as the investment threshold). Using (|2.t)|l and 
H3.17|l . we find for x < h: 



V+ix;h) = {1 - q)iU%u){x) 



/3+(A+-/3+) 
A+ 

X+ 



e-^^y[{l - q)-'n-{l)qM{l)Ge-+y - C]dy 

x+y>h 

K-{l)qM{l)G 



1 



+ 00 



A+ 



h — x 

e''K-{l)qM{l)G _ C_ 
(l-<j)(/3+-l) "/?+ 



dy 



h—x 



Ce'^ ydy 



(3.19) 



3.3 Timing an investment 

In the case of the stream g{x) ~ e^qM {1)G — {1 — q)C, H2.11|l assumes the form 
(1 - g)-ie^K-(l)gM(l)G - C = 0, (3.20) 



and the investment threshold is 



(l-g)g 

Kqil)qM{l)G' 



(3.21) 



After h* is found, the EPV of the project can be calculated with the help of 
((2.6(1 . In particular, if the probability density is given by ((3.8(1 . we can use ((3.19(1 
with h — h* , and substituting 1(3.21(1 , find 



(\+ _ g+)„P+{x~h') 



(3.22) 
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3.4 Timing an exit 



Consider a special case when the profit function in Subsection 2.4.2 is given as 
f{Xt) = Ge^* . Then g{x) is a linear combination of exponentials, so that we 
can calculate U'^g{x) using (|3.13|) : 

Uy{x) = q^C - (1 - q)-'K+{l)qM{l)Ge-, 
and find the optimal disinvestment threshold from l|2.15|l : 

._^il-q)C 



K+{l)qM{l)G 



(3.23) 



After h^, is found from H3.23|l . we can calculate the expected present value of the 
gains from scrapping the production unit at time 0, V~ {x; h^), using H2.8|l and 
ijT^ . Similarly to lE23l, for x > K, 



V^{x;K) ^ Cq 



Finally, the value of the firm is the sum of the EPV of the perpetual stream of 
profits plus the option value of scrapping: for a; > /i,, 

pxn (\- _ f>-\ (3- {x-h,) 



4 Incremental capital expansion 

In this Section, we assume that the production function depends only on capital, 
and that G{K) is differentiable, concave, and satisfies the Inada conditions. 
A similar situation was considered in [1] for a two-period model of partially 
reversible investment, in |14| for the geometric Brownian motion model, and in 
0] for Levy processes. At each time period t, the firm receives e-^*G{Kt) from 
the sales of its product, and, should it decide to increase the capital stock, suffers 
the installation cost C ■ [Kt+i — Kt)- The firm's objective is to chose the optimal 
investment strategy K. = {Kt+i{Kt^ Xt)}t>ii Kq = K, Xq = x, which maximizes 
the NPV of the firm: 

ViK,x) = sup^;-^ 

Here we treat the current log price x and capital stock K as state variables, 
and AI^ as a sequence of control variables. Due to irreversibility of investment, 
Kt+, > Kt, Vt. 

To ensure that firm's value H4.1|l is bounded, we impose a resource constraint: 
there exists K < oo, such that Kt < K, \ft. The resource constraint, condition 
1)3. l|l . and properties of the production function ensure that the value function 
H4.1|l is well defined. 



J2q'ie'''G{Kt)-CiK. 



t+i 



Kt)) 



(4.1) 
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Formally, the manager has to choose both the timing and the size of the 
capital expansion. However, it is well-known (see, for example, ^3]) that for 
each level of the capital stock, it is only necessary to decide when to invest. 
The manager's problem is equivalent to finding the boundary (the investment 
threshold), h{K;C), between two regions in the state variable space {K,x): 
inaction and action ones. For all pairs {K, x) belonging to the inaction region, it 
is optimal to keep the capital stock unchanged. In the action region, investment 
becomes optimal. To derive the equation for the investment boundary, suppose 
first that every new investment can be made in chunks of capital, AK, only''. In 
this case, the firm has to suffer the cost CAK, and the EPV of the revenue gain 
due to this investment can be represented in the form of the EPV of the stream 
giXt) = qM{l){G{K + AK) - G{K))e^' - (1 - q)CAK. On the strength of the 
result of Subsection 3.3, the optimal exercise boundary is determined from the 
equation U'^g{h) = 0, which can be written as 

(1 - qy\M{l){G{K + AK) - G{K))K-{l)e'' = GAK. (4.2) 

Dividing by AK in (|4.2|l and passing to the limit, we obtain the equation for 
the optimal threshold, h* — h*{K): 

qM{\)K-{l)G'{K)e^ = C(l - q). (4.3) 

Equivalently, the optimal exercise price is 

qMil)^q{l)G'iKy 

The rigorous justification of this limiting argument can be made exactly as in 
the continuous time model in 4 . Let h — h{K;A) be the solution to H4.2|l . 
Then the option value associated with the increase of the capital by AK, at the 
price level e^, is 

(1 - q) (71 1 (x) (^^M^{G{K + AK) - G{K))k- il)e^ - CAK^ . 

As AK ^ 0, we have h = h{K;A) ^ h*{K); therefore, dividing by AK and 
passing to the limit, we obtain the formula for the derivative of the option value 
of future investment opportunities w.r.t. K: 

V°^\K,x) = (l-g)t/|^l[,.,+^)(a;) (l)e- - . (4.5) 

Substituting G from 14.4|l into (|4.5I) and using the definition oi k^{1), we obtain 



V°^\K,x) = il-q)E 



J2Q'G'iK)e^'\Xo^O 



t=o 



*The authors are indebted for this simplifying trick to Mike Harrison; our initial proof was 
more involved. 
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The last formula factors out the contributions of the infimum and supremum 
price processes to the marginal option value of capital. The first expectation on 
the RHS decreases if the probability of downward jumps in prices increases, and 
the second expectation increases if the probability of positive jumps in prices 
increases. Hence the marginal option value of capital increases in downward un- 
certainty and decreases in upward uncertainty. The overall effect of uncertainty 
is ambiguous. 

Now the marginal, or shadow, value of capital is given by 

Equation (|4.t)|) expresses the marginal value of capital as the difference of two 
components. The first one is the expected present value of the marginal returns 
to capital, given that the capital stock remains constant at the level K in the 
future. The second component is the marginal option value of the future invest- 
ment opportunities. This value is subtracted because investing extinguishes the 
option. A similar result was obtained in ,1' for a two-period model of partially 
reversible investment. 

Exactly the same calculations which lead to 13.22|l allow us to derive from 

vriK..)^cL_^4^_. („) 

Substituting (|4.4|1 into the last equation, we derive 
Integrating (|4.8|l w.r.t. capital, we obtain 



The value of the firm is V{K, x) = e^G(A')/(l - q_M{l)) + V°p\K, x). Suppose 
that the capital stock available for investment, K, is very large, and function 
G'{KY^ is integrable on [l,-|-oo). The n we can obtain a simpler formula by 
replacing the upper limit K with -|-oo. In the case of Cobb-Douglas production 
function G{K) = dK^ , we have G'{K) = ddK^~^, therefore the integrals con- 
verge iff {9 — l)^'^ < —1, or, equivalently, /3+ > 1/(1 — 9). If this condition is 
satisfied, the value of the firm for x < h*{K) is 

ViK,x) ^ 



qM{l) 

X+ -13+ fn-{l)qM{l)\^* ^i-/3+g;3+.^i-/5+(i-e) 



(/3+-l)A+V (l-q) J {d9r+{p+{l-9)-l) 
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One of the standing assumptions in corporate finance is that the value of a firm 
follows a (geometric) Brownian motion. Our solution clearly shows that the 
firm's value is not described even by a Markovian process. 

Notice that the proof of (|4.3|) in 0] was based on the reduction to the case of 
the perpetual American call, and therefore the generalization for more general 
dependence on the stochastic factor was not possible. Here the result holds for 
any continuous increasing revenue fiow R{K, x), and the formula for the optimal 
investment threshold obtains in the form: 

UyTRK{K,h)^C, (4.9) 

where the EPV-operator [/^ and transition operator T act w.r.t. the second 
argument (we have qTRk instead of Rk because the revenues will start to accrue 
the period after the investment is made). For instance, if the firm faces the 
operational cost a + bKe^'^^^, then the revenue fiow is R{Kt,Xt) = e-^*G{Kt) — 
a — bKe'^'^^^, and instead of 14.4(1 . we now have 

(1 - q)-'[qM{l)^-{l)G'{K)e'^ - bqM{l/2)^; {l/2)e'^/'] - C = 0. (4.10) 

The function on the LHS in H4.10|l changes sign only once, and therefore the 
solution to equation H4.10|l gives the optimal investment threshold. 



A Proof of optimality of h* and 

On the strength of Lemma on p. 1364 in to prove the optimality of the 
solution V{h*;x), it suffices to check the following two conditions: 

V{h*;x)>m&x{G{x),0}, V x, (A.l) 

and 

{I -qT)V{h*;x)>0, Vs. (A.2) 

By our choice of h* , 'w{x) > for all x > h* , therefore V{h*;x) = (1 — 
q)U'^l[h'^+oci)w{x) > for all x. Using 1(2. 6|l . we represent V{h*;x) in the 
form 

Vih*;x) = {l^q)U'^Uy{x)-{l-q){U'^l^_^^^,)w)ix) 
— Uxgix) — v{h*;x) = G{x) — v{h*;x), 

where v{h*;x) = (1 — q){U'^l(^-oc,h')w){x). Due to the choice of h* , we have 
w{x) < Vx < ft-*, so that v{h*;x) < for all x, and hence V{h*;x) > G{x) for 
all X. We conclude that ((A.ip holds. 

Under a very weak regularity assumption on p, it is proved in [7] that for any 

h, 

V{h- x) = qTV{h; x), x < h. (A.3) 

Using more elaborate arguments as in [H] and [7] in continuous time, one can 
prove that ((A. 3(1 holds if the probability density p exists. Thus, ((A.2(l holds 
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on {—oo,h*), and it remains to verify ljA.2|) for x > h* . Introduce W{x) = 
{I ~ qT)V{h*;x). Using ijOJ and lETj), we obtain 

W{x) = {I~qT){l-q)U^^l[h'^+oo)Ul9{x) 

= ({4)-4[,,,+^)C/l5(x), (A.4) 

and also 

W{x) = {I-qT){l~q)Um,',+o.)Ul9(x) 

= (/ - qT){l - q)U%Uj,g{x) -{I- qT)il - q)U'll^_^^^,^U'^g{x) 



g{x) + {-I + qT){l - q)Um^^,n^)Uj,g{x) 



But if a function u vanishes outside (—00, /i), then 



!7|w(x) = E 



^q^u{x + Xt) 



.t=o 



0, M x>h, 



as weh. Therefore, for x > h* , 



Wix) = 5(2;) + 9(1 - q)TU]^li-oo.h*)Uy{x) 



(A.5) 



By our standing assumption, g is non-decreasing, hence from (|A.5|I . W{x) is 
non-decreasing on [h*,+oo). We apply [/^ to l|A.4|l and obtain 



J^q'Wia 



(A.6) 



Suppose that W{h*) < 0. Then there exists hi > h* such that W{x) < for all 
X G {h*,hi). It follows that for the same x, the RHS in (|A.6p is non-positive. 
But for these x, the LHS is positive by the very definition of h* . Hence, our 
assumption W{h*) < is false, and since W is non-decreasing on [h* , +00), the 
condition HA.2() follows, and the proof of optimality is finished. 

The verification of the optimality conditions ljA.l|) - (|A.2p for the exit problem 
is quite similar to the proof above: just replace X,2(^, C/^, and [h*,+oo) with 
X, X,f7^,C/^ and (— cx),/i,], respectively. 



B Transition densities given by exponential 
polynomials 

B.l The case of three exponentials 

In this Subsection, we demonstrate how to obtain the transition density of a 
desired shape. The density 13.8|l has a kink (and maximum) at the origin. If we 
want to have a smooth p (and allow for the maximum to be not at the origin), 
we need to use more than two exponential functions. Suppose that we want to 
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model a density which has the maximum on the positive half-axis. Then we use 
one exponential on the negative half-axis, and two on the positive one: 

(c+A+e-^^- - c+X+e-^inMo,+oo){x) ~ e'^' -l^^^,o]{x), (B.l) 

where c^,c^ and c~ are positive, and < < A]^ < Aj. Later in this 
Subsection, we show that for any choice of A^ < < Aj^ < A^, equation ljB.l|l 
with c^, cf , C2 given by simple formulas (jB.6p defines a probability density, 
which has the maximum on the positive half-axis. See Figure 1 for an example. 
Similarly, one can construct a 3-parameter family of probability densities which 
have the maximum on the negative half-axis. Should one wish to have a smooth 
probability density which has the maximum at the origin, one must use two 
exponential functions on each half-line or exponential polynomials of the form 
{ax + b)e^^. 

The moment generating function of the probability density (|B.1|) is 

C+A+ C+A+ c-A- 
M{z) - ^ ^ 22 



Aj-z AJ-z A - z 

At the end of this Subsection we will show that the fundamental rational function 
1 — qM{z) = has 3 roots, all of which real. Call these roots /3~, (3^ and /?^. 
We have 

A- < /J- < < /?+ < A+ < A+ < /?+. (B.2) 
Clearly, 1/(1 — qM{z)) can be represented in the form 

1 '^t ^2 fT3 n\ 

1 - qM{z) = Jl^z ^ ^ 'P^~z' ^ 

where a+j = l/((?Af'(/3+)), j = 1,2, and a' = l/{qM'{l3~)), and therefore, 

poo 

U],g{x) = a+ e-^^yg{x + y)dy (B.4) 



e ^ 



00 

PtVnf^ 4- .Aril, - n- I f^-f^^Vr 



e ' ■^yg{x + y)dy - a I e ' yg{x + y)dy. 
<y —oQ 

Similarly, we can consider a probability density given by linear combinations of 
two or more exponents on each of the half-axis. If we use two exponentials for 
each, we have two roots pf , j — 1, 2, of the "characteristic equation" l—qM{z) = 
on each half-axis, and l|B.2p - (|B.4|l change in the straightforward manner. One 
can also use more than two exponentials on each axis, and obtain more elaborate 
probability densities. 

Now we show that any choice A^ < < Aj^ < Aj defines a probability 
density. Three conditions: p{x)dx = 1, p is continuous at 0, and p is 

smooth at 0, give a linear system of three equations 

cf-c+ + cr = 1, 

c^Xt -C+X+ +C-X- = 0, (B.5) 

ct{xtr-4i^tr+c-{x-r - o. 
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Using Cramer's rule, it is easy to find a unique solution {cf , cj, c ) to ljB.5|l for 
any A~ < < Aj: 



~A-A+ 



(B.6) 



(A+-A+)(A+-A-)' 
+ ^ ^A-A+ 

(A+-A+)(A+-A-)' 

(A+-A-)(A+-A-)- 

It is easily seen that cf , 62 and c~ are positive, and that p is positive as well. 

The roots of 1 — qM{z) are found as follows. Clearly, 1 — qM{z) has 3 roots 
at most. As z ^ A^ + 0, 1 — qM{z) —00, and the same holds as z ^ Aj^ — 0, 
and as z ^ A;^ + 0. Under condition g e (0, 1), 1 — qM{0) — 1 — q > 0, and 
1 — gM(+oo) = 1 > 0. Hence, on each of the intervals (A",0) {0,\i ), and 
(A^, +00), 1 — qM{z) changes sign. Therefore, on each of these three intervals, 
there is exactly one root, which we have called and respectively. 

B.2 General scheme for the computation of [/^ and 

Step 1. Calculate the moment-generating function M(z), and consider the 
rational function 1 — qM(z). Find the roots of the denominator, A^, and the 
numerator, , with their multiplicities (sign "+" for the roots on the positive 
axis, sign "-" for the ones on the negative axis). 
Step 2. Define 

j 3 k ' 

step 3. If all the roots (3^ are simple, we represent and Kq{z) in the 

form 

<(z) = 4(oo)+^-^, ^-(,) = «-(oo)-5]-^, (B.9) 
where 

The case of multiple roots can be treated similarly. 

Step 4. For a continuous g satisfying H3.5|l . we can calculate 



r+00 ^ 

Uj^gix) ^ K+{oo)g{x) + y2(^k / e'^i'^gix + y)dy, 

I. Jo 



(B.IO) 
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and for a continuous g satisfying we can find 

U^gix) = K- {^)g{x) +Y.^k e-^" ^9{x + y)dy. (B.U) 

k 

Under condition H3.7|l . all the roots [3^ are outside [(T~,(t+], therefore H3.5f) and 
H3.5|l ensure the convergence of the integrals in (|B.10|I and (jB.ll|l . 



B.3 Value of investment project 

By assuming that the probability density is given by exponential polynomials 
on each half-axis and using (|2.6|) and (|B.10|I . we find for x < h: 

V+{x;h) = 

-K-{l)qM{l)G 



k 



^-0^ y+x+y^y _ (J / g-/3, y^y 



h — x 



h—x 



e''K-{l)qM{\)G _ C_ 



If the threshold is chosen optimally, then we use H3.21|l and find 



(B.12) 



(B.13) 



B.4 Value of scrapping 

After /i, is found from H3.23|l . we can calculate the expected present value of the 
gains from scrapping of the firm at time 0, F~(a;;/i*), using p. 8(1 and H3.23|l . 



Similarly to l|B.13|l . for a; > ft.*, 

Finally, the value of the firm is the sum of the EPV of the perpetual stream of 
profits plus the option value of scrapping: for x > h^, 



V{x) 



i^)^ '\i^-Pk)Pk■ 



l-qM{l) 



B.5 Capital expansion 

Exactly the same calculations which lead to H3.22|l allow us to derive from 1(4. 7|l 



y < ^^Uc'iKft (B 14) 
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Integrating (|B.14|I w.r.t. K, we obtain 



Y (/?; - \ {^-i)c ) Jk ^ ' 

Suppose that the available capital stock, is very large, and the functions 
G'{K)^k are integrable on [1, +00). Then we can obtain a simpler formula by 
replacing the upper limit K with +00. In the case of Cobb-Douglas production 
function G{K) = dK^ , we have G'{K) = d6K^~^, therefore the integrals con- 
verge iff {9 — < —1 for all k. Let be the smallest positive root, then the 
equivalent condition is > 1/(1 — 9). If this condition is satisfied, the option 
value of investment opportunities is 



iPt - V (1 - 9) ; {d9)f'i (/?+(! -9)-l) 



C Expected waiting time 

Assume that the spot log-price x is less than h*, and consider the waiting time 
till the investment is made. If the transition density is given by H3.8|) . then 
the expected waiting time is finite iff A+ + < 0, or, equivalently, 

TO = = 1/A- + 1/A+ > 0, (C.l) 

and if (|C.ip holds, then 

E[R.,] = -{h* -X + 1/X+). (C.2) 

TO 

Condition ljC.l|) has a clear interpretation: the expected waiting time is finite 
iff the drift of the log-price, to, is positive, and if it is positive, then (|C.2|I 
says that the expected waiting time is inversely proportional to the drift. It is 
also proportional to the distance to the barrier plus the positive term which is 
independent of the distance. 



To obtain the formula l|C.2|) for the expected waiting time, we need to cal- 
culate the limit in H2.17|l . Using H3.17|l . we obtain for x < h*: 

R+(\+ _ R+\ /• + °° 

Up[H',+oo){x) = (l-q)-' ' ' ^ e-P^y^j,,,+^){x + v)dy 

^1!%^'^^^ = (1 - q)-'e-P'('^''^\l - /?+/A+). 
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Hence, we have to calculate the limit 



lim 



1 ~e^' 



■(3+{q)ih'--x) e-/3+(9)(''*-^)/3+(q)/A+ 



1-q 



Both terms are positive, and if /3"'"(g)/(l — q) is unbounded as g ^ 1, the limit 
is clearly infinite. From (|3.1U|I . we find that this possibility does not realize iff 
A" + A+ < 0. If + A+ < 0, we obtain 



and (|C.2|I follows. 
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Figure 1: Density IjB.ip . Parameters: = — 5,A]^ = 5,A^ = 7.5 
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